Superlattice formed by quantum-dot sheets: density of states and IR absorption 
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Low-energy continuous states of electron in heterosrtucture with periodically placed quantum-dot 
sheets are studied theoretically. The Green's function of electron is governed by the Dyson equation 
with the self-energy function which is determined the boundary conditions at quantum-dot sheets 
with weak damping in low-energy region. The parameters of superlattice formed by quantum-dot 
sheets are determined using of the short-range model of quantum dot. The density of states and 
spectral dependencies of the anisotropic absorption coefficient under mid-IR transitions from doped 
quantum dots into miniband states of superlattice strongly depend on dot concentration and on 
period of sheets. These dependencies can be used for characterization of the multi-layer structure 
and they determine parameters of different optoelectronic devices exploiting vertical transport of 
carriers through quantum-dot sheets. 

PACS numbers: 73.21. Cd, 73.21.La, 78.67.Pt 



I. INTRODUCTION 

Heterostructures formed by quantum dot (QD) sheets 
are widely investigated and used in different devices, such 
as lasers, photodetectors, and solar cells, see for re- 
view. In such heterostructures, not only the additional 
localized states of electrons captured into QDs should 
be taken into account but also the continuous electronic 
states, which are subjected to reflections on periodically 
placed QD sheets, should be modified significantly. Such 
a periodical perturbation gives rise to a superlattice (SL) 
with energy spectrum formed by gaps between allowed 
minibands. In contrast to the standard case, Q an 
additional damping of electronic states takes place due 
to scattering on inhomogeneties of QD sheets stemming 
from a random in-plane distribution of QDs. But such a 
damping appears to be weak for low-energy region. As 
a result, SL effect should be essential near the edge of 
interband absorption in host material, which is propor- 
tional to the density of states of SL, or under IR transi- 
tions from doped QDs into miniband states. To the best 
of our knowledge, these phenomena were not considered 
based on a simultaneous description of SL minibands and 
damping effects in spite of the structures under consid- 
eration are routinely used in different optoelectronic de- 
vices. At the same time the opposite case of 3D ordering 
of the closely spaced QDs, when SL is formed as a result 
of the tunneling mix between intra-QD states, was ana- 
lyzed [f| and demonstrated experimentally, see @, 0] and 
references therein. Because of this, it is important and 
timely to develop an adequate theory of low-energy elec- 
trons interacting with the periodically-placed QD sheets 
and to study the optical response of SL which can be used 
for characterization of structures under consideration and 
for description of different optoelectronic devices. 

In this paper we study low-energy electronic states, 
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with energies in the vicinity of the conduction band ex- 
tremum, in heterosrtuctures formed by QD sheets of pe- 
riod I using the effective-mass equations for the Green's 
function averaged over randomly placed QDs in each 
sheet. In contrast to the standard theoretical descrip- 
tion based on the averaging over 3D or 2D space, Q, 
here we perform the averaging over QD sheets with the 
identical statistical characteristics. As a result, we ob- 
tain the inhomogeneous along SL axis Dyson equation 
where the self-energy function can be replaced by the 
boundary conditions at QD sheets. Since the damping of 
the low-energy states is weak, one can consider SL which 
characteristics are determined by an effective potential 
localized at the sheet positions, z = nl, n = 0, ±1, .... 
The strength of this potential is determined by the con- 
centration of QDs and shape of QD potential. With re- 
spect to low-energy states, QD can be considered as a 
short-range defect (which are widely investigated during 
past 50 years, see Refs. 9) if the low-energy interval un- 
der consideration is smaller than the QD binding energy. 

The density of states in SL depends on the period I and 
on the parameter determined by a strength of QD's po- 
tential described within the short-range approximation. 
Spectral dependencies of interband absorption between 
the heavy-hole and SL states are proportional to the den- 
sity of states in c-band. In addition, the anisotropic ab- 
sorption coefficient, originated due to mid-IR transitions 
from the doped QD ground state into the miniband states 
of SL, is obtained through the QD concentration and the 
SL parameters. We found that the efficiency of mid-IR 
photoexcitation is comparable to the contribution of wet- 
ting layers formed under QD sheets [l(| if doping levels 
are the same. But the spectral dependencies are very 
different for these two mechanisms. Thus, it is demon- 
strated that the results obtained can be used for charac- 
terization of structure under consideration. It is more im- 
portant that the SL parameters determine a mechanism 
of vertical transport for underbarrier electrons which is 
a key process in different optoelectronic devices exploit- 
ing multi-QD sheets. Similar mechanism of transport 
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through vmderbarrier states of IR photodetectors formed 
by GaAs/AlGaAs-based SL was considered in Ref. 11. 

The paper is organized as follows. In Sect. II we de- 
scribe the model of periodical sheets formed by randomly 
placed QDs and evaluate the Green's function averaged 
over random positions of QDs. SL effects on the density 
of states and on the process of anisotropic photoexcita- 
tion of QDs are considered in Sect. III. List of assump- 
tions used and concluding remarks are presented in the 
last section. Appendix contains the justification of the 
effective SL approach employed in the calculations per- 
formed. 



II. MODEL 

The electronic states near c-band extremum of het- 
erostructure, which is formed by QD sheets placed in 
host material, are described by the effective mass Hamil- 
tonian 
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where p is the 3D momentum operator, m is the effective 
mass, and u(r— R rJ t) is the potential energy of QD placed 
at coordinates R r fe = (x r k,rl). Here r labels sheet (r — 
0, ±1, ±2, . . .) placed with the period I and k stands for 
position of QD over the rth sheet given by 2D random 
coordinate x r fc (k — 1,2 ... N where N is number of QDs 
over each sheet with the normalization area L 2 ). Electron 
of energy E is described by the Green's function Qe (r, r') 
governed by the equation 



E + iX-H)Q E (r, r') = 6 (r - r' 



(2) 



with A — » +0 and the 3D 5- function S(Ar). Below we 
consider the averaged over all QD positions Green's func- 
tion Gfi (r,r') = (Qe (r, r')) where the averaging over rth 
sheet is performed according to [8j 
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and (...) includes the averaging over all sheets. 

Using the (p, ^-representation (p is 2D momentum) 
one obtains the Dyson equation governing the averaged 
Green's function as follows 

Ge p (z,z') — gE P {z - z') (4) 
+ dzi dz 2 gE P (z - zi^Ep (zi, z 2 ) G Ep (z 2 , z') ■ 



Here gEp (z — z') is the free Green's function which is 
governed by Eq. (2) with the Hamiltonian p 2 /2m, so 
that 



g Ep (Az) = - 
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if e p > E and the imaginary factor iwE — e v should be 
used in (5) if e v < E. Within the self-consistent Born 
approximation, the self-energy function He p (z\, z 2 ) in 
Eq.(4) is given by 
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where u(q, z) is the 2D Fourier transform of u(r) and 
tiqd is the QD concentration over sheet which does not 
dependent on r, i.e. we consider identical QD sheets. 

Further, we restrict ourselves by the low-energy region 
where scattering on a QD can be described by the short- 
range potential u(r) ps UA(r) with the form- factor A(r) 
localized in volume ~ a 3 (a stands for the characteristic 
size of QD). We also neglect high-order corrections to the 
self-energy function (6), see below the diagram expansion 
of Fig. 4 and discussion in Appendix. Since the kernel (6) 
is located near QD sheets with zi,2 ~ rl and the Green's 
functions vary over scales h/y/2m\E — e p \, the integral 
equation (4) is transformed into the finite-difference one: 



G Ep {z,z') = g Ep {z - z') 
-A Ep 2J 9e p (z - rl)G Ep (rl, z'). 
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The self-energy function (6) is written here through the 
factor 
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which is the same for any QD sheet [we moved S r . . . from 
Eq. (6)toEq. (7)]. Instead of Eq. (7), one can determine 
Ge p (z, z') from Eq. (2) with the free Hamiltonian p 2 /2m 
and describe the QD sheet effect adding the boundary 
conditions 
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at sheet positions z — rl. This result was evaluated after 
acting of the operator E + iX — p 2 /2m on the integral 
Dyson equation (4) and subsequent integration of the 
intergo-differential equation obtained over the QD posi- 
tions (rl - 0,rl + 0). 

Within the second-order Born approximation we use 
GE P (rl,rl) ~ gE P (0) in the self-consistent equation (8), 
see Ref. 8 for detais, and the momentum-independent 
factor As in Eq. (9) takes the form 



A E = A [l + i 
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Here we estimate A for the case of short-range defect 
within the Kostcr-Slater approach [l2[ and ~p E is the 
3D density of states which is taken at the cut-off en- 
ergy e a ~ (irh/ 'a) 2 /2m. Since E <C e , damping of low- 
energy states is weak and one can replace the complex 
boundary condition (9) by the effective potential energy 
— A^2 r 5 a (z — rl) with 5 a (A.z) localized in the interval 
\Az\ < a, so that in the framework of the effective SL 
approach Ge p (z,z') is governed by the one-dimensional 
equation: 

(E + i\ - e p - H ± ^j G Ep (z, z') = S(z - z'), (If) 

i 

2m 



with the electron effective mass in the GaAs matrix, 
to. Thus, the Green's function is expressed using the 
standard relation Q between Ge p (z,z') and the solu- 



tions of the eigenstate problem for SL, [13j H±iJj 
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The last equation determines the dispersion 



relations s np± and the cigcnfunctions i]j 
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Here p± is 



quasimomentum (\p±\ < nh/l), n labels minibands, and 
the wavefunction takes form 



4^ = VV, 
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where the reflection coefficient and the normalization fac- 
tor, R np± and ip np± , are expressed through p± and k np± . 
14] The energy s np± — (Tik np± ) 2 /2m is founded from the 
dispersion equation 



cos — — = cos k 
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which is written through the characteristic wave vector, 
K = Am/h 2 ~ n 3 nQ£,a/2, see Ref. 13 for details. 

The dispersion relations for lower minibands deter- 
mined by Eq. (13) are shown in Fig. 1 for dimensionless 
parameter Kl varied between 1 and 8 when the trans- 
formation from the weakly-coupled SL (if Kl <2) to 
the tight-binding regime of coupling (if Kl >4) takes 
place. The characteristic energy e; = (irU/l) 2 /2m is 
about 3.2 meV for SL of period I =40 nm. For SL 
formed by InAs QDs embedded by GaAs matrix Kl «3.1 
if tiqd — 5 x 10 10 cm~ 2 . As a result, minigaps exceed 
5 meV for the tight-binding regime, see Figs, lc and Id 
when dispersion laws are close to cosine and sine depen- 
dencies, for odd and even n respectively. For the weakly- 
coupled SLs the dispersion laws are formed by parabolic 
curves modified near p±l/Ti = 0,7r with gaps ~1 meV, 
see Figs, la and lb. In contrast to SL corresponding to 
under-barrier tunneling regime, }13| if Kl <1.5 one ob- 
tains the lowest miniband at finite p±l/h only, as it is 
shown in Fig. la. This is because of absence of solution 
for Eq. (13) at p± — > and k np± l <C 1. Such a peculiar- 
ity change the density of states and the edge of mid-IR 
absorption if Kl < 1.5, see Figs. 2a and 3a below. 
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FIG. 1: Miniband energy spectra, E/ei versus p±l/h, of the 
effective SL determined by Eq. (13) for Kl =1 (a), 2 (b), 4 
(c), and 8 (d). 



III. RESULTS 

Using the model described above, we consider in this 
section the density of states in SL formed by QD sheets, 
and calculate the absorption coefficient under mid-IR 
photoexcitation from ground levels of doped QDs into 
miniband states of SL. 



A. Density of states 

The density of states is introduced through the aver- 
aged Green's function by the standard formula [8] 



Pe = ~^ Im J dr (Ge (r, r)) 



(14) 
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where 2 is due to spin degeneracy and L 3 is the normal- 
ization volume. The lower expression is obtained for the 
case of negligible damping in Eq. (10) using of the ef- 
fective SL approach determined by Eqs. (11)-(13), see 
the energy spectra plotted in Fig. 1. The integration 
of (^-function over p gives the 2D density of states, p 2 D, 
and after integration of ^-function over p± the density of 
states should be replaced by constant if E belongs to nth 
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FIG. 2: Normalized density of states PeI/P2d versus E/ei 
given by Eq. (15) for the parameters used in panels (a-d) of 
Fig. 1. Dotted curve in upper panel corresponds to the 3D 
density of states oc \f~E if SL effect is negligible, A — > 0. 



gap. In nth miniband (below nth gap) , the integral over 
p± should be taken over the interval (0,pe) where pe is 
found as a root of the equation E — Sn PE . As a result, 
Pe takes the form: 



PE 



Pip 
I 



1 



n, 

-PeI/(ttH), 



E C nth gap 
E C nth band 
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and a shape of pe is determined by the gap- induced steps 
with transitions between them determined by the mini- 
band dispersion laws. 

In Fig. 2 we plot the dimensionless density of states, in 
units P2d/1, for the same parameters as in Fig. 1. For the 
weak coupling regime, the jump of pe at E — > appears 
due to the cut-off of the lowest miniband at finite p±l/h, 
c. f. Figs, la and 2a at E/ei < 2. With increasing of Kl 
under transition to the tight-binding regime, the energy- 
independent gap contributions to the density of states 
increase and pe between these steps is transformed from 
oc \J~E dependency shown by dotted curve in Fig. 2a to 
the arccosine dependencies. In addition, the bottom of 
lowest subband is shifted to energies ~ ei. Since pe is 
connected directly to the shape of interband optical spec- 
tra, see Ref. 13, the step-like dependencies permit one 
to extract Kl value which determine the bandstructure 
of SL according to Eq. (13). 

Let us compare the energy scale of SL effect, deter- 
mined by and the SL effect due to the wetting layer 
contribution analyzed in Refs. 10. For the parameters 



given at the end of Sect. II, one obtains that the contri- 
bution of QD sheet with tiqd = 5 x 10 10 cm~ 2 is reduced 
~2 times in comparison with the wetting layer effect if 
levels of electron doping are the same. Thus, an interplay 
of both mechanisms should take place for tiqe, > 10 11 
cm~ 2 . For such a case, the interband optical spectra 
should be dependent on both the QDs contributions and 
the wetting layer contributions. 



B. Photoionization 



The anisotropic absorption coefficients aJi and arc 
determined from the general Kubo formula as follows: 



y/ecujL 6 4-' 

do' 



x |(£|e||,_L • v\6')\ S (e s ~ e s > + huj) , (16) 
where e is the dielectric permittivity of the host semi- 

I 1 2 

conductor and the matrix element (£|eii j_ ■ vja) cor- 
responds to transitions between 5- and ^'-states of en- 
ergies sg and Eg' under radiation with polarization orts 
eiM- We use the equilibrium distributions f(eg) — » 1 
and f(es + fvui) because the only localized states are 
populated at temperatures lower the binding energy \Eq\. 
Due to the in-plane isotropy of the problem, we separate 
the cases of s- and p-polarized radiation corresponding 
to the polarization orts and e z . Neglecting the over- 
lap between QD states and taking the ground state wave 
functions 'F p in the momentum representation (P is 3D 
momentum) we transform Eq. (16) into: 



xe 



J 



i(Pr- P'r')/ft 



4ne 2 



^p*p' / dr f dr' 
pp' •* 



(e||P)(e||P') 

PzPt 



K A 



(r',r). (17) 



The contribution of miniband states are described here 
through the average of the exact Green's function 
QE(r',r) with the exponential factor corresponding to 
random QD positions (here Ap = P — P'): 



^A P , £ (r',r) = (^ e i;A P R -/ ri Im£ B (r',r) 



(18) 



which is analyzed in the Appendix. Within the low-order 
approach, the correlation function (18) takes the form: 

K Ap Ar', r) « N QD jS Apfl lmG E (r', r), (19) 

where NqoL/1 is the total number of QDs in the nor- 
malization volume I? and the averaged Green's function 
Gs(r',r) was considered in Sect. II. 

Using the ground state wave function typ written in the 
Koster- Slater approach [l2[ and neglecting the damping 
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hAco/s l 

FIG. 3: (Color online) Spectral dependencies of dimensionless 
absorption coefficients determined by Eqs. (20) and (21) for 
the same conditions as in Figs. 1 and 2. Solid and dashed 
curves correspond to the perpendicular and parallel polar- 
izations, respectively. Dotted curves correspond to the case 
A — > 0, when SL effect is negligible. 



correction in Eq. (10) we transform Eq. (17) as follows 



(2ire) 2 riQ D 
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dze -zp ± z/h^(np ± ) 



S (hAuj 



Here P = (p,p_i_) and we have replaced Ge p (z', z) from 
Eq. (19) using the wave function (12). In the expres- 
sions for a\i' ± integrals over p-plane and over z are taken 
analytically and the spectral dependencies of IR absorp- 
tion are obtained after the double numerical integrations 
over the transverse momenta p±_ and pj_ . The dimension- 
less spectral dependencies are plotted in Fig. 3 for the 
same conditions as in Figs. 1 and 2. The characteristic 
absorption ao is given by 



a 



{Ae) 2 n QD I ,-, 



(21) 



c^/em\E \/2 V^o. 
and ag ^3 cm for the above listed parameters. Thus, 



for the maximal absorption, when TiAu>/ei ^20 - 30 
or HAuj ~ l-^o |) one obtains a max ^45 cm" . Since 
ali'' ± cx hqd/I 4 , the maximal absorption increases up 
to a max > 10 3 cm -1 if tiqd > 10 11 cm -2 and I ~20 nm; 
an approximation of low QD concentration remains valid 
for such a set of parameters. Further increase of a max 
is possible in the case of heavily doped SL, with a few 
electrons captured in QD. 

Anisotropy of absorption is about 20% without any 
strong dependency on effective potential, c. f. Figs. 3a 
- 3d where parameter Kl varies from 1 to 8. Peculiari- 
ties of miniband spectra are visible clearly in a^j starting 
from Kl >2 while at does not show any peculiarities 
at the edges of minibands. This is due to different se- 
lection rules for transverse and longitudinally polarized 
excitations: in the last case, transitions are forbidden at 
edges of minibands and the spectral dependencies rmain 
smooth. In addition, Fig. la shows a jump of a^ at 
%Auj = which is similar to the jump of the density of 
states in Fig. 2a (we do not consider IR transitions into 
shallow underbarrier states at TiAuj < 0). In Figs. 3b-d, 
shifts of absorption edges to finite TiAlj > take place 
due to lower miniband shifts, see Figs, lb-d and 2b-d. 



IV. CONCLUSIONS 

In summary, we have developed the theory of the su- 
perlattice formed by periodically placed quantum dot 
sheets. It was found that the damping due to random 
in-plane positions of dots is weak and effect of the sheets 
on electronic states can be described using of the effec- 
tive boundary conditions. Within this approach we have 
demonstrated that the miniband density of states, which 
describes the interband absorption, and spectra of mid- 
IR photoexcitation of doped quantum dots into mini- 
bands strongly depend on parameters of quantum dot 
sheets. Visible anisotropy of the absorption coefficient is 
also found, with transverse absorption which is strongly 
modulated by the miniband spectrum of SL. 

Now we discuss the main assumptions in the calcula- 
tions performed. We restricted ourselves by the vicin- 
ity of c-band using the effective-mass approach in Eq. 
(1) and in further consideration of the photoionization 
process. In order to describe the energy intervals com- 
parable to the gap, one needs to use the multi-band 
kp-Hamiltonian for more detailed description of QD 
states. [l5| We consider the case of low QD concentra- 
tion (uqd/1 ~ 10 15 cm~ 3 in our numerical estimates) 
and the electron-electron interaction effect on the energy 
spectrum; thus, the IR-absorption should be weak. Nu- 
merical estimates for the SL parameters were performed 
here based on simplified description of QD as an isotropic 
short-range defect with the binding energy correspond- 
ing to typical QD. This approach gives approximate SL 
parameters only and a more precise description should 
be based on a numerical solution of the self-consistent 
Dyson equation taking into account a real potential of 
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P P ■ P /^~~ \ P 

z z' z z' z z^__^/z 2 z' 

FIG. 4: Self-consistent Dyson equation for averaged Green's 
function Ge p (z,z') and the self-energy function Esp (zi, z%) 
shown in upper and lower lines, respectively. 

QD. [l|, [l5| Because parameters of QD sheet (materials, 
concentration, and shape of QD) can be very different, 
such a consideration should be performed for different 
specific cases (e.g. for Ge/Si-based or Aj/Byj-based QD 
sheets, for review see Ref. 16). 

To conclude, we believe that the results obtained will 
stimulate an investigation of underbarrier vertical trans- 
port of carriers in order to verify SL effect on electronic 
properties of structures formed by QD sheets. The spec- 
tral and polarization dependencies of the mid-IR pho- 
toexcitation are convenient for direct measurements be- 
cause the valence band states are not essential. These 
results should be important for description of different 
devices utilizing periodical QD sheet structures. 
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Appendix 

In order to estimate the corrections beyond the effec- 
tive potential approach used in Eqs. (9)-(ll) we consider 
here the method of calculations in more details. Us- 
ing the pz-representation, one obtains the self-consistent 
Dyson equation (4) for the averaged Green's function 
Ge p (z,z') shown by a bold line as it is plotted in Fig. 
4. Within the second-order Born approximation, we use 
the free Green's function in the self-energy function (6) 
given by the first diagram of the set for E^ p shown in 
the lower line of Fig. 4. The next corrections in this set 
can be neglected under the standard condition [8J 

£»|£ Bp |~A (A.l) 



and we arrive to Eq. (7) using the free Green's function 
in E^p determined by Eq. (6). 

More complicate consideration is necessary for the cor- 
relation function Ka p e(y, r') appearing in Eq. (17) be- 
cause of the random factor exp(iApR rfc /?i) describing 
positions of QDs. Instead Eq. (18) it is convenient to 
consider the generalized expression 




FIG. 5: Diagram expansion for correlation function 
ICa p e(v, r') written through the diagram set for vertex part 
shown in lower line. 

which is shown in Fig. 5. Here a dotted curve corre- 
sponds to the averaged factor 

( cxpf--^Ap-R rifcl ju(r-R r2fc2 ) V (A.3) 

\r 1 k 1 r 2 k2 I 

while dashed curves in Figs. 4 and 5 stand for the 
paired QD potentials. After summation over all reducible 
diagrams, Kape(*, r') is written through the averaged 
Green's function and the vertex part which is given by 
the set shown in the lower line of Fig. 5 with the initial 
vortex determined from (A.3) as follows 

7A P (ri, r 2 ) = n QD ^ u (~~7~> Zl ~ rl ^j ( A - 4 ) 

xe -|(Apx 1+rI)i !)^ ri _ r ^ 

The first correction to Eq. (19) appears, if we use (A.4), 
as the vertex part in the diagram expansion for correla- 
tion function shown in Fig. 5. Performing straightfor- 
ward calculations under the condition (A.l), one obtains 
that this correction and next contributions are negligible 
in comparison with Eq. (19). 
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